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We study thermalization by applying gradient expansion to the Kadanoff-Baym 
equations of the 2PI effective action to two-loop in a theory with Dirac fermions 
coupled to scalars. In addition to those chemical potentials which equilibrate in the 
on-shell limit, we identify modes which are conserved in this approximation, but 
which relax when off-shell effects are taken into account. This implies that chemical 
equilibration does not require higher loop contributions to the effective action and is 
compatible with the gradient expansion. We explicitly calculate the damping time- 
scales of both, on- and off-shell, chemical equilibration rates. It is shown that off-shell 
equilibration is suppressed by the thermal width of the particles in the plasma, which 
explains the separation of on- and off-shell chemical equilibration time-scales. 



I. INTRODUCTION 

Kadanoff Baym equations are an effective method for describing the out-of-equilibrium 
dynamics of quantum fields. For many particular models and initial conditions, thermaliza- 

ly from numerical solutions to these equations, see e.g. 



tion has been demonstrated successful 
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refs. [I|,l2|,l3|, y.ay.0 , 18|, 19|, HO, 
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13| for some works in the field. While details can de- 
pend on the particular settings, it has been revealed that starting from far-from-equilibrium 
initial conditions, thermalization proceeds typically in several usually overlapping stages: 
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• Dephasing/prethermalization an In this stage, a constant equation of state is 
reached and an approximate equipartition between kinetic and potential energy is 
attained. This dephasing effect does not rely on scattering-driven processes and is a 
pure quantum phenomenon. 

• Kinetic equilibration , Q: Within this stage, the quasi-particle distribution func- 
tions approach the Bose-Einstein or Fermi-Dirac form, respectively. The interactions 
which drive the system towards kinetic equilibrium do not necessarily relax chemi- 
cal potentials through violation of particle number. For example, this is the case for 
elastic two-by-two scatterings. 

• Chemical equilibration j^, Q : Finally, processes that do not conserve particle numbers 
eventually relax the chemical potentials. 

This last stage of chemical equilibration is the subject of the present study. The inclu- 
sion of off-shell effects constitutes one qualitative difference between the Kadanoff-Baym 
and the Boltzmann equations. This has been mentioned as a possible explanation for the 
discrepancy between numerical solutions within the two approaches j^]. Besides, Kadanoff- 
Baym equations include memory integrals and treat deviations from equilibrium to all order 
in derivatives, usually referred to as gradients. Boltzmann equations do not incorporate 
memory integrals and are first order in gradients. 

For systems which are close to equilibrium, gradient expansion and neglect of memory 
integrals is justifiable. On the other hand, due to the possible presence of quantities which 
only equilibrate when off-shell effects are included, it is desirable to account for these within 
an analytical approach. For the case of scalar X(j)^ theory in 2 -|- 1 dimensions, Boltzmann 
equations have been generalized to include off-shell processes. Numerical solutions to these 



equations are found to be in accordance with the solutions to the fu 



tions, but in disagreement with the on-shell Boltzmann equations 



1 Kadanoff-Baym equa- 



i5i,y,H]. 



The explanation 



is that while 2^2 scatterings lead to kinetic equilibration at a time-scale proportional to 
A^, chemical equilibration is driven _by 1 3 off-shell processes that are suppressed by 
higher powers in A; see also e.g. ref. 9] for related discussions and refs. j^, [l^ for a numer- 
ical analysis in 3+1 dimensions. The case of an interacting system of fermions and scalars 



is numerically treated in 12|, where solutions to the Kadanoff-Baym and the Boltzmann 



equations are compared. It is found that there are chemical potentials which are conserved 
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by the Boltzmann equations but which are relaxed when Kadanoff-Baym equations are ap- 
phed. This phenomenon is explained by the fact that the quasi-particle approximation of 
Boltzmann equations leads to spuriously conserved quantities. 

The aim of the current work is to show how to reproduce on- and off-shell chemical 
equilibration in a semi-classical approach, without resorting to numerical solutions. As a 
concrete model, we consider a Dirac fermion coupled via a Yukawa coupling to a complex 
scalar. We expect that this model shares the essential features for the discussion of chemical 



n 



equilibration in the quark-meson model considered in refs. j^, [1^ and the models usually 



considered for leptogenesis [16| and for chargino mediated electroweak baryogenesis [17| • The 
classical Boltzmann equations can be derived from the Kadanoff-Baym equations by using 
several approximations. These are the gradient expansion (to first order), the expansion 
in coupling constants (typically up to the first nontrivial order that accounts for scattering 
or decay processes) and the quasi-particle (or on-shell) approximation. We will show that 
by omitting the last approximation, full chemical equilibration can be described by semi- 
classical transport equations. In particular, the gradient expansion is valid in this regime 
and we present analytic expressions for the on- and off-shell relaxation time-scales of the 
chemical potentials. 

The paper is organized as follows. In the next section we set up notation and in section UTTl 
we introduce the gradient expansion. In the two subsequent sections, the model and the 
Boltzmann equations in the on-shell limit are discussed. Finally, we calculate the various 
time-scales of chemical equilibration in section IVIt before we conclude in section IVIII 



II. THE KADANOFF-BAYM EQUATIONS 



In this section, we derive the Kadanoff-Baym equations and set up notation. We follow the 
nomenclature of ref . 18| , where a more complete treatment can be found. Statistical systems 
in equilibrium and non-equilibrium can be described by the Kadanoff-Baym equations which 
are the Schwinger-Dyson equations in the in- in formalism and in Wigner space The 
in-in-formalism is designed to describe the time-evolution of density matrix elements rather 
than of scattering processes to which the usual in-out formalism applies. The in-in formalism 



is implemented by integrating along the Schwinger-Keldysh 20|, \2m\ closed time path (CTP) 
from a finite time t = to to t = +00 and back. The Green functions for bosons and fermions 
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are then defined as the path ordered operators 

iA{u,v) = {n\Tc<j){u)<f)\v)\n), (1) 
iSai3{u,v) = {Q\TciJa{u)'4!p{v)\Q) . (2) 



Depending on the locations of the coordinates u and v on the branches of the CTP, we 
obtain the four Green functions 

iA (m, v) = {n\T(f){u)4>^{v)\n) , (3a) 

iA<{u,v) = {n\(f)^{v)(f){u)\n) , (3b) 

iA>{u,v) = {n\(f){u)(j)\v)\n) , (3c) 

iA\u,v) = {n\f(j){u)(j)'^{v)\n), (3d) 

and similar expressions for the fermions. In the following, we focus on the bosonic functions. 
The fermionic functions fulfill similar corresponding equations (see ref. 18|). From the 
definitions ([3]), it is obvious that only two of the four functions are independent while the 
other two can be expressed as 

A\u,v) = e{uo-vo)A>{u,v) + e{vo-uo)A<{u,v), (4a) 
A\u,v) = e{uo-vo)A<{u,v) + e{vo-uo)A>{u,v). (4b) 

Besides, the Green functions have the following hermiticity property 

(iA<'>(M,t;))^ = iA<'>(^;,u). (5) 

The retarded and advanced propagators are defined as 

A'' = A* - A< = A> - A*, (6a) 
A" = A* - A> = A< - A*, (6b) 

which can be used to define the spectral and hermitian functions 

^0 = ^(A^-A'^) = ^(A>-A<), (7) 

Ah = ^{A'- + A'') = -i sign(Mo - vo)A^. (8) 

The usefulness of these definitions becomes obvious for the thermal equilibrium Green func- 
tions in Wigner space, where one obtains for a free boson of mass M 

A^{p,x) = j d^pe'P 'A^ix + r/2,x- r/2) = 7r5{e - M^) sign(fco), (9) 
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and 

= 2rf^{k,) A^, = 2{rQ{ko) + I) A^, rf^{k,) = (10) 

Using this notation, the Kadanoff-Baym equations for the bosons in Wigner space read [1^ 

{nl ± ir^) * A^''^ = 1, (11) 

* A<'> - n<'> * A;, = ^(n> * A< - n< * a>), (12) 

where 

A* B = Ae-^'^^^^"-^'^^^ B (13) 

denotes the Moyal star product and the various functions 11 are self-energies which have the 
same time-ordering prescriptions as defined for the Green functions in eqs. ([3]). The self- 
energies can be deduced from the two-particle-irreducible (2PI) effective action. Besides, we 
have defined 

nl = k'-M'-u,, r^ = ^(n>-n<). (i4) 

In the following, we discuss this system of equations in the gradient expansion. 



III. GRADIENT EXPANSION 

We consider spatially homogeneous and isotropic systems, such that all functions depend 
on the time t = xq, but not on the spatial coordinates x. At rather late times, when 
the system is not too far from equilibrium, one expects that the relaxation of the system 
towards equilibrium is driven by interactions which allow for a gradient expansion in the 
weak coupling regime.^ Thus, one can expand the Moyal star product f|T3l) 

A*B^AB- ^{A, B},, + 0{dl), (15) 

where we have defined the Poisson brackets {A, B}ph = d^AdkB — dkAd^B. 
At first order in the gradient expansion, eq. (fTTj) reads 

{nl ± ir,^) A^'" - '-{nl ± ir^, a^-'^},, = i, (le) 



This argument can however be jeopardized by oscillation effects in the case of several flavors as discussed 
in ref. 
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which is solved by 

Here, we have used the fact that {A, 1/A}pb = —{A, A}pf,/A'^ = 0. This leads to the spectral 
and hermitian functions 

r fl^ 

The remaining information is encoded in the anti-hermitian part of eq. f[T^ . which reads 

- {Q\ iA<^>U + {in<'>, iAJp, = n> A< - n< A>. (19) 
Since we assume that the coupling constant y is small, the self-energies are at late times 



effectively of first order in the gradient expansion. We therefore neg. 



ect the self-energies in 



the Poisson brackets (in contrast to the methods used e.g. in ref. [15|]), which leads to 

2A;o<9tiA<'> = n>A< - n<A>. (20) 
Analogously, one finds for a free fermion in equilibrium 

iS< = -2n'^^{ko)A^, iS> = 2{l-n'^^{ko))A^, n'^{ko) = (21) 

G ~t~ J- 

and the fermionic Kadanoff-Baym equations 

{f}^ ± iU) ^u, ^^'^^5 = 1' (22) 

7odtiS<'> = ^(^>5< - p<S>) + h.c, (23) 

where ^ denotes the fermionic self-energy and fl^ = j) — and = ^(^^ — ^^). This 
implies that 

where we have defined = 2^^ ■ = {^_4, ^^}, = ^^^^ and fi^ = fi^fi^. In 
general, the fermionic pro pagator has a complicated structure and e.g. contains a pole from 
a collective excitation [241]. However, the corresponding residuum is exponentially small for 
most particles in the plasma such that one can safely ignore this subtlety. 
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IV. THE MODEL AND ITS CONSERVED CHARGES 

For our purpose of an analytic calculation of the separation of time-scales of on- and 
off-shell chemical equilibration, we consider a theory with a Lagrange density as in ref. 13] 

L = iipp^ + df,(j)^d^(j) - (j)^M'^(f) -y^{PR(f) + PL(j)^)tlj, (26) 

with a small real coupling constant y. Using the two-loop 2PI effective action, this La- 
grangian leads to the following expressions for the self-energies [l^ 

iU<'>{k,x) = -y'l ^^S{k-p-q)TT[Pr,iS>'<{-p,x)PLiS<'>{q,x)], (27) 

and for the fermions (^ = + ^i?) 

-^Sik -p- g)iA>'<(-A;, x)P^i5<'>(-g, x)Pn, (28a) 

-j^Sik -p- q)iA<'>{k, x)PniS<'>{-p, x)Pl. (28b) 

We now define the charges 

Qs = J -0^,2koiiA< + 2ei-ko)A^), (29a) 

= J -0^^2k'^koiiA< + 2ei-ko)A^), (29b) 

^Tr7oPL/ii(i5< + 2e{-qo)A^), (29c) 

Pl/r = - J ^jTq^loPL/R{iS< + 2e{-q,)A^), (29d) 

where the terms involving the spectral functions A have vanishing time-derivatives and have 
been introduced to make the expressions finite. The interpretation of these quantities in the 
on-shell limit is given in the next section. When using the kinetic equations fl20p and (1231) . 
the form of the collision terms immediately leads to the following conserved quantities: 

Qr + Ql = const. (30a) 
2Qs-Qr + Ql = const. (30b) 
Ps+Pr + Pl = const. (30c) 

We emphasize that these conservation laws are due to the symmetries of the collision terms 
and make no assumption about the Green functions and their spectral properties. Of course, 
this is reflecting that above charges are conserved due to Noether symmetries with respect 
to rephasings and translations of the Lagrangian (|26|) . 
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V. THE ON-SHELL LIMIT 

The on-shell limit arises when neglecting the influence of interactions on the spectral 
functions. Notice, that in order to obtain the correct results, one has to take into account 
that the Breit-Wigner width behaves as F — > sign(fco)e for both bosons and fermions. In 
this case, we recover the spectral functions of the free theory 

A^ik.x) 7r5(A;2-M^)sign(A;o), (31) 
A^ik,x) ^ 7r^5(A:2)sign(fco), (32) 



that fulfill 



^^OoU A 1 

— 2koA^ = 1, 

ZTC 

^270^^ = 1. (33) 



Next, we employ the ansatz 



iA< = 2nsiko)A^, (34a) 

iA> = 2ins{ko) + l)A^, (34b) 

iS< = -2{nLiko)PLA^PR + nn{ko)PRA^PL), (34c) 

iS> = 2iA^ - nL{ko)PLA^PR - nR{ko)PRA^PL) , (34d) 

with the generalized particle distribution functions n(fco). This turns the conserved quanti- 
ties in eq. (l29i) into 

/d^k 

j^{ns{uj)-nsi^)), (35) 

/d^k 

^^{u;ns{u;)+u;ns{u;)), (36) 

/d^k 

j^^{nL,R{k) -nL,R{k)), (37) 

PIr = I 0^,{kn,,R{k) + kn,,R{k)), (38) 

where we have introduced the bosonic dispersion relation cj^ = + M^. Besides, we have 
defined the anti-particle distribution functions for negative energies 

ns{uj) = -{ns{-uj) + 1), nifiiko) = 1 - nL,R{-ko). (39) 
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From these definitions, it is clear that the conserved quantities represent total energy 
and charges. In particular, the knowledge of the initial charges and total energy should 
suffice to determine the final equilibrium state of the system that is characterized by the 
temperature and the chemical potentials after equilibration. Assuming kinetic equilibrium, 
the Kubo-Martin-Schwinger (KMS) relation including chemical potentials are 

A> = e^(P»-'^«)A<, PlS> = Sp'-^'^^PlS^, PrS> = -e^(f«-^«)PR^<. (40) 

These relations then imply that the distribution functions are of the form 

n^_^ = l/(e^(po-'^^.«) + l), ns = l/{e'^^"'^ - I). (41) 

Under the close-to-equilibrium assumption that fiL,R,s/T ^ 1, we can linearly relate the 
charges to the chemical potentials by introducing statistical factors as 

Ql,R,S = -^^L ji sfJ-LAS ! (42) 
D 

where 

oo 

kL,R=l2T-^ I - nf) = ^ / -'d-J-^. = ^ ' (43) 



oo , 

'^=^'^"/ (^""^-B^+ = i / ^^^r' - j^j^^^- 

M/T 

In the limit M/T 0, the latter integral evaluates to ks = 2, and we plot the statistical 
factors function of M/T in FigE 

Next, we note that the self-energies fl27l) and fl28l) inherit the KMS property fl40l) from 
the Wightman functions, 

U>{k) = n<(fc)e^('=«"^^+^'«), (45a) 

= -mk)S''-''^''\ (45b) 
f>L{k) = -^<(A;)e^('=«-^«-^s). (45c) 

From the Feynman rules in the Schwinger-Keldysh formalism (in particular from energy- 
momentum conservation) it is possible to derive that these relations even hold to arbitrary 
order in perturbation theory. Consequently, we find for the collision terms 

U>{k)A<{k) - U<{k)A>{k) = U>{k)A<{k) (1 - e^('^^-'^«-'^s)) ^ (45^) 

mk)sm - mk)snk) = (i - e^^^--^^^^^^) , (46b) 
mk)sm - r^msm = mk)sm (i - e^^^-^-^^^o , (46c) 
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such that these terms vanish and a static solution of the transport equations in eqs. (120|) 
and fl23|) is obtained if 

/is + A^R - yUL = 0. (47) 

Hence, one would expect that the final state is determined by the temperature and the three 
conserved quantities in eq. (l30ll only, and that it satisfies the equilibrium condition (HTll 
above. 

However, it turns out that in the on-shell approximation, not only the combinations fl30|) 
are conserved, but also additional charges in the particle and anti-particle sectors sepa- 



rately 12|. This is a purely kinematic effect. Consider a scalar decaying into a fermion/anti- 
fermion pair. If all three particles are on-shell, all energies have to have the same sign. This 
is most easily seen in the rest-frame of the scalar, where ko = ±M. In this frame, the 
momenta of fermions differ only by a sign and hence their energies are equal and therefore 
% = Po = ^o/2. This means that the collision integrals in eqs. (l20l) and (l23l) have in the 
on-shell limit only a support in the regions 

sign/co = signpo = signgg- (48) 

This is sufficient to show that the number of fermions plus twice the number of scalars is 
separately conserved for particles and for anti-particles, such that the Boltzmann equation 
cannot reproduce chemical equilibration in the on-shell limit. Formally, we can express these 
conservation laws by defining 

/d'^k 
^{nL,R,s{^) + nL,R,s{^)) (49) 

and noting 

^Qs + Ql + Qr = const. (50a) 
Ql ~ Qr = const. (50b) 

In order to distinguish from the conserved charges defined in eqs. ( !30ll . we refer to these 
combinations as total particle number densities. We introduce an additional set of chemical 
potentials Jil,r,s to account for total particle numbers different from their equilibrium values. 
These are defined through the relations 
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PlS> = -e^iPo-pL^isnpo)p^g<^ p^s> = -e^(Po-fin^isnpo)p^g<_ (52) 

Now by taking combinations of ^l,r,s and piL,R,S) "we can adjust the number of particles and 
anti-particles independently. 

For a support of the form as in eq. (H8l) . the collision terms for non- vanishing fiL,R,s read 

n>(A;)A<(A;) - n<(A;)A>(A;) = n>(A;)A<(A;) (l - e^'S'^('=o)'3(^^+'^«-^^^)) ^ ^53^^^ 

mk)sm - mk)sm = mk)s^ik) (i - e^^^-^'^^w~-,---^^+-^^^) , (53b) 

mk)Sm - mk)Sm = mk)S^ik) (1 - e^^^-^'^^^(--^^--^-^-^^^) , (53c) 
and these vanish as long as 

fis - ^L- P'R = 0. (54) 

Notice that this is only true in the two-loop approximation of the effective action. 

To summarize, we have defined six chemical charges, Ql,r,s and Ql,r,s- Through on-shell 
processes, two linear combinations of these charges are not conserved and equilibrate accord- 
ing to the conditions (H71) and (|5^ . Therefore, in the on-shell limit there remain four con- 
served chemical charges as defined by eqs. f l30al I30b"l) and ( l50i) . While the charges f l30al I30b"|l 
are conserved by virtue of Lagrangian symmetries, the conservation of the charges fl50|) is an 
artifact of the on-shell approximation. Notice that this argument is in principle not based 
on the gradient expansion, since for these time-independent chemical potentials the collision 
terms vanish to all orders in gradients as long as the support of the collision integral is of 
the form as in eq. f HHl) . 

VI. TIME-SCALES 

The considerations of the last section show that the relaxation of some chemical potentials 
is described by on-shell Boltzmann equations. The rate of these relaxation processes is 
given by the thermally averaged matrix elements of the on-shell scattering process. Off-shell 
equilibration on the other hand is additionally suppressed by the Breit-Wigner width of the 
particles in the plasma, which is why we expect these processes to be slow compared to 
the on-shell relaxation. In the following, we analytically calculate the on-shell relaxation 
time-scale Ton and the off-shell time-scale scale Toff in linear response theory, and we show 
that indeed, they exhibit a separation Ton ^ Tog. For this purpose, we assume that the 
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deviations from equilibrium are of the form 

6A<{k) = 6A>{k) = 2(36fis{k)n^{k){n''^{k) + 1)A^, (55) 
6S<{k) = 6S>{k) = -2(3{PL6fiL{k) + PR6fiR{k))n'^''{k){l - n'^{k))A^. (56) 

These expressions follow from the definitions flMl) with 

= e/3(fco-5Ms(fc)) _ 1 ' ^lA^) = g/3(fco-5A.i,fl(fc)) + 1 ' (^^) 

when expanded for 5ixs,r,l/T <^ 1. 



A. Charge Equilibration 

Let us now calculate the relaxation rate for the charge densities that is damped by on-shell 
processes. Therefore, we consider the case of chemical potentials without energy dependence, 

Using the kinetic equations (l20i) and (l23l) . one obtains in linear response three equations 

dtQL = ^C°"(/i5 + /iii-/iL), (58a) 
dtQR = -^C°'^(/^5 + fiR- fii), (58b) 
dtQs = -^C-^'if^s + fiR- /iL), (58c) 



what is consistent with the conserved combinations in eq. (1301) . The equation for the non- 
conserved combination of chemical potentials (1471) is 

dtif^s + m- I^l) = + ^) C°°(/i5 + m- I^l) , (59) 



where in on-shell approximation 



= -y^l ^^SlJls{k-p-q)TT[PmS>{-p)PLiS<{q)]zA>{k)] (60) 

X sign(go)5(g^) sign(po)<5(p^) sign{qo + po)5{{q + pf - M^) . (61) 
After energy and angular integration, this reads 

= y'M' J '^,e{Aqp - M')nf{p)nf{q){l + rf^{p + q)) . (62) 



13 




0.01 



M/T 



FIG. 1: The plots show the coefficients /cl,_r and % as functions of the scalar mass M. 



The symmetry of this expression becomes more apparent when noting that 



nf{,p)nf{q){l + n^{p + q)) = (1 - nf{p)){l - nf{q))n^{p + q) 



(63) 



The integral (162|) has an interpretation in terms of processes that describe the decay of the 
scalar into a fermion-anti-fermion pair and annihilation of a fermion-anti-fermion pair into 
a scalar. These processes are suppressed for small scalar masses, since this reduces the phase 
space of the kinetically allowed final states. We plot C°° in Fig. [21 from where its features 
can be verified. 

From there the typical time-scale of charge equilibration can be read off: 

6 



7 



Q 



The resulting time scale for the choices y'^ = 1.0 and = 0.1 is shown in Fig. [31 
For M/T ^ 1, it is instructive to evaluate the relevant integrals as 



(64) 



, — _M l'M\ 2 

ks ~ ^yl8/n^e t | _ 



and 



such that 
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FIG. 2: The plot shows the collision terms C°° and as functions of the scalar mass M. The 
combination jy^ has a small residual dependence on y and we use y = 1.0 in the plot. 
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FIG. 3: The plots show the time-scales of on-shell chemical equilibration, 7*^ given by eq. ()64p . 
and the three different time-scales given by the eigenvalues of the system in eq. (j69p as a function 
of the scalar mass M. In the left (right) plot a coupling y^ = 1.0 (y^ = 0.1) has been used. 

Therefore, chemical equilibration according to (H7j) is attained faster for larger M. We 
emphasize however, that this does not imply that physical interaction rates become fast. The 
opposite is true, as can be seen from the exponential decrease of C""^. Chemical equilibration 
can occur, because in the large mass limit, a change in the chemical potential /x^ only 
corresponds to an exponentially small chan ge i n the physical charge density Qs- The increase 
of 7*^ with larger M has been noted in refs. 25|,l26|]. There it has been shown that this feature 



15 



simplifies considerations of cfiemical equilibration for electroweak baryogenesis and for the 
conversion of baryon-minus-lepton number to baryon number, as relevant for leptogenesis. 

B. Total Particle Number Equilibration 

In order to investigate the equilibration of total particle number, consider chemical po- 
tentials of the form 

SfiL,R,s{k) = fiL,R,ssign{ko). (68) 

Using this ansatz in the kinetic equations (l20l) and ( |23l) and linearizing in the chemical 
potentials then leads to 

'^ksdtfis = - 2fi+) - r^fis, (69a) 



^kL,Rdtfi+ = C°"(/i5-2/i+), (69b) 
6 

^ kL,Rdt^- = -C°^Ji-. (69c) 
6 

where we have defined the combinations /i± = |(/iL ± JIr) and the off-shell decay rate 

= d'pd^qd^k ^^^ _ ^ _ g)Tr[P^^g>(-p)P^zg<(g)]iA>(A:) 

x[sign(A;o) - (sign(go) + sign(po))/2]^ (70) 

Notice that without the off-shell contributions, there would be two unsuppressed modes, 
namely /2_ and ksfis + ki^jifl^. On the other hand, including the off-shell effects and in the 
limit of /C°^ <^ 1, these are still approximate eigenmodes which are damped at the rates 

7f-|J^ZiC°^ 7? - (71) 

In addition, there is the approximate mode fls — 2/i+, which is suppressed by on-shell 
equilibration at the rate 

l^^^,{2kl^ + k-s')C^\ (72) 

Let us now proceed with the evaluation of the decay rate close to equilibrium. The 
main contributions to this integral come from the regions, where two of the particles are 
on-shell, but the third is off-shell violating the condition in eq. ( HHi) . We hence replace two 
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of the Breit-Wigner spectral functions by delta functions. For simplicity, we present only 
the contribution where the scalar is off-shell, in which case one obtains the contribution 

roS _ „,2 /" d^vd^q 



/d pd Q 

X sign(go)(5(g^) sign(po)5(p' 



(p2 + 2p ■ g + g2 _ M2)2 + r^po + go, \p + gl) 
X [sign(po + go) - (sign(go) + sign(po))/2]l (73) 

Notice that if the two fermions are on-shell with sign(po) 7^ sign(go), the momentum of the 
scalar = {u, k) is necessarily space- like, u"^ — k"^ < 0. The corresponding Breit-Wigner 
width is close to equilibrium given by 

-{k+uj)/2 



T^{uj,k) = y^{uj^-k'^) / ^n^{p) 

J{k-c,)/2 Stt/c 



Stt 

Integration over the energies yields 




(74) 



rod 2 I dpdq sin 9d9 ^ ^ 

= y / pg(l-cos^) 

X (%''(-p)%'^(g)^B (P - g) - n^{p)np{-q)nT^i(l " P)) 
r<^(b - g|, y/p2 + g2 _ 2pqcos6) 



(75) 



(2j9g(l - COS^) + M2)2 + r2(|p - g|, ^p2 + g2_2pgcOS^) ' 

Here, we have used that go and po need to have opposite signs. Notice that this does not 
contain the potentially harmful pole in the bosonic distribution function, since the width F 
vanishes in this case. A plot of C^^ is provided in Fig. [2l 

Since the off-shell effects are proportional to the Breit-Wigner width, they are propor- 
tional to y'^ and hence suppressed compared to on-shell effects for small couplings. The 
resulting equilibration rates given by 7*^ as in eq. and by the eigenvalues of the system 
in eq. (169!) for the choices y^ = 1.0 and y^ = 0.1 are shown in Fig. [31 

It is again possible to make an estimate for the M/T ^ 1 limit. The integrand in eq. fl75|) 
receives only relevant contributions from p,q ^ M, such that we can replace the nominator 
by M^. Numerical evaluation then yields 

~ 0.126 1/^^. (76) 
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The suppression for large M is not exponential, but given by the Breit-Wigner width. From 
Fig. [3] it is however evident that when M/T is not too large and there is a sizable number 
of scalar particles S present in the heat bath, the on-shell equilibration rate is dominating 
over the off-shell processes. Since the on-shell rate is oc and the off-shell rate oc y^, this 
separation of scales becomes more pronounced for smaller y. 

VII. DISCUSSION AND CONCLUSIONS 

We have calculated the various time-scales of chemical equilibration in a theory with a 
Dirac fermion coupled to a scalar. For that purpose, we have applied a gradient expansion 
and a linear response ansatz to the Kadanoff-Baym equations. The different time-scales 
arise hereby from processes that are driven either by on-shell or off-shell effects. While on- 
shell equilibration is described in a quasi-particle picture with 5-functions accounting for the 
on-shell conditions, we have shown that by including finite width effects, the calculation can 
straightforwardly be generalized to determine the rate of off-shell relaxation. We emphasize 
that already at leading loop order in the 2PI effective action, a self-consistent solution to 
the spectral and kinetic equations leads to off-shell chemical equilibration. In addition, the 
gradient expansion is consistent with off-shell chemical equilibration. Parametrically, we 
find for the time scales of on- and off-shell chemical equilibration 

7on OC ?/^ 7ofr oc ^/^ (77) 

what explains the separation of time scales. In our model, equilibration of charge results 
from on-shell processes, while there exist two linear combinations of total particle numbers 
that only relax through off-shell contributions. 

For our discussion of the off-shell relaxation rate, we have focused on the contributions to 
the collision integral where the scalar particle is off-shell. Contributions of similar size (and 
same sign, as can easily be checked) are expected from the domains of integration where 
the fermions are off-shell. The numerical values presented should hence not be taken at 
face value, but might be larger by a factor 2 — 3 in the full calculation, when taking the 
fermionic off-shell effects into account. Since we find that the rate of off-shell effects is of 
order y^, three-loop contributions to the effective action might also be equally important 
as the off-shell effects arising from the two-loop effective action. However, the analyses of 
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Kadanoff-Baym equations in the literature we referred to also take only the two-loop order 
of the effective action into account. Hence, we neglected higher loop orders in the effective 
action to facilitate a direct comparison. 

We would like to stress that our analysis does not only confirm the common lore that 
off-shell effects can be important in certain models. In the model at hand, full equilibration 
requires the violation of eq. (HHI) . which is a stronger condition than deviation from the on- 
shell limit. As a consequence of eq. fHHj) in conjunction with the Breit-Wigner form of the 
spectral function, we see from eq. (176|) the decaying behavior of the off-shell equilibration 
rate for large M. In contrast, if the spectral function would be Gaussian with a width 
proportional to y^, off-shell equilibration would be exponentially suppressed. 

The advantage of transport equations using the gradient expansion in Wigner space is 
that they do not contain memory integrals and therefore are much easier to solve numerically 
than the full Kadanoff-Baym equations in coordinate space even if no further approximations 
as e.g. the on-shell approximation or an ansatz in terms of chemical potentials are used. 
Hence, it is important to understand the range of applicability of transport equations. In 
summary, our results indicate that the use of the gradient expansion is justified to describe 
thermalization at rather late times (namely after effects as prethermalization and damping 
are concluded) if the off-shell effects are properly taken into account and coupling constants 
are rather small. This is seen explicitly in the presented analytical approach: In the on- 
shell approximation, the collision terms vanish at all orders in gradients in some cases even 
though the system is not in its true equilibrium state, while once including off-shell effects, 
the system achieves full equilibration already at first order in the gradient expansion. Hence, 
higher orders in the gradient expansion are not decisive to describe full equilibration at late 
times qualitatively in the present system. 
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